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I analyze the algebraic patterns underlying the structure of scattering amplitudes in quantum 
field theory. I focus on the decomposition of amplitudes in terms of independent functions 
and the systems of differential equations the latter obey. In particular, I discuss the key role 
played by unitarity for the decomposition in terms of master integrals, by means of generalized 
cuts and integrand reduction, as well as for solving the corresponding differential equations, 
by means of Magnus exponential series. 


1 Introduction 

High energy particle collisions are the ideal framework for accessing new informations on matter 
constituents and forces of nature. The higher the energy of the colliding particles, the richer the 
landscape of the produced ones. The discovery of new physics interactions cannot be disentan¬ 
gled from the discovery of massive, heavy particles, emerging from collisions of ever increasing 
energy. On the other side, by increasing energy, also the probability of producing many light par¬ 
ticles is enhanced. Therefore, advances in High Energy Particle Physics necessarily depend on 
our ability to describe the scattering processes involving many light and heavy particles at very 
high accuracy, hence they depend on our capability of evaluating Feynman diagrams. Beyond 
leading order (LO), Feynman diagrams represent challenging multidimensional/multivariate in¬ 
tegrals, whose direct evaluation is often prohibitive, therefore the computation of scattering 
amplitudes beyond the LO is addressed in two stages: i) the decomposition in terms of a basis 
of functions, and ii) their evaluation of the elements of such a basis, called master integrals 
(Mis). In this contribution, I elaborate on the algebraic properties of Feynman integrals, which 
can be exploited for decomposing them in terms of Mis and for computing the latter. The 
techniques I discuss can be applied to generic amplitudes, and have a impact on high-accuracy 
prediction for collider physics, as well as for exploring the more formal aspect of quantum field 
theory. 

Let us observe that amplitudes can be decomposed in terms of independent functions, ex¬ 
actly like a vector can be decomposed along basic directions. One needs a basis and a projection 


technique. The latter is necessary to extract the coefficients of the linear combination. Fac¬ 
torization is the basic idea we are going to elaborate on. Factorization is ubiquitous in the 
discovery of new mathematical and physical concepts. Complex numbers emerged from factor¬ 
izing the simplest number we may think of, i.e. 1 = (—i)i m , quantum mechanics relies on the 
factorization of the identity matrix, I = |n)(n|; Dirac equation emerged from factorizing the 

d’Alambertian operator, i.e. □ = (—ip)(ip). What does happen when amplitudes factorize? 

Cutting a virtual particle and bringing it on the mass shell ( p 2 = m 2 ), turned out to be a suit¬ 
able projection technique yielding amplitudes decomposition. Why multiple-cuts are important? 
First, because multiple-cuts yield functions identification. Since any diagram is characterized 
by its internal lines, a given master diagram is univocally identified by a cut-diagram where 
all internal particles are on-shell. Moreover, when applied to amplitudes, multiple-cuts behave 
like high-pass filters, which isolate only the diagrams that have those internal lines to be cut, 
while the others are automatically discarded. Therefore, by considering all possible cuts of an 
amplitude, in a top-down procedure, from the maximum number of cuts to the lowest one, it is 
possible to build a (triangular) system of equations from which all coefficients can be determined. 


2 Integrand decomposition 



Figure 1 - Tree-level recurrence re¬ 
lation. 


Tree-level scattering amplitudes obey a quadratic recurrence 
relation 1 (BCFW), depicted in Fig.l, whose derivation relies 
on Cauchy’s residue theorem. Since tree amplitudes are ratio¬ 
nal functions of kinematic variables, the BCFW recurrence can 
be understood as due simply to partial fractioning 2 , because 
residue theorem applied to rational functions amounts to partial 
fractions. Is that just accidental, and holding for tree-level am¬ 
plitudes, or partial fractioning can be exploited also at higher 


orders? 

The integrand reduction algorithm 3 had a dramatic impact on our ability of computing 
one-loop amplitudes. The basic idea lies in the exis¬ 
tence of a relation between numerators and denom¬ 
inators of scattering amplitudes which can be used 
to decompose the integrands of one-loop amplitudes 
in terms of integrands of Mis. The amplitude de¬ 
composition in terms of Mis is then achieved after 
integrating the integrand decomposition. The co¬ 
efficients of the Mis are a subset of the coefficients appearing in the decomposition of the 
integrands. Therefore, within the integrand reduction algorithm, coefficients can be determined 
simply by algebraic manipulation, with the great advantage of bypassing any integration. 

The idea behind the GoSam framework 4,5 is to combine 
automated diagram generation and algebraic manipulation 
with the integrand-level reduction, implemented in Samu¬ 
rai 6,7 and Ninja 8,9 and the tensor decomposition of Golem95 
The code is very flexible and it has been employed in several 
applications at NLO QCD accuracy, studies of BSM scenar¬ 
ios, electroweak calculations, and recently also within NNLO 
calculations. It is interfaced to several MonteCarlo event 
generators, like Sherpa, Herwig, aMC@NLO. GoSam was 
used to evaluate the NLO QCD correction to pp -A Hjj, Hjjj 
(in the infinite top-mass limit) 11,12 , w hich required an ex¬ 
tension of the integrand decomposition methods 8 . The evaluation of the virtual amplitudes 
for pp —> Hjjj has been further optimized, enhancing the numerical accuracy and reducing 
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Figure 3 - LO and NLO total cross sec¬ 
tions for pp —> #+1,2,3,jets at LHC (13 
TeV). 




Figure 2 Multiloop integrand decomposition. 







the computing time 13 , ending up into a new phenomenological analysis 14,15 , obtained with the 
tandem of Go Sam and Sherpa, see Fig. 3. 

The extension of the integrand decomposition beyond one-loop has been proposed in 16 , and 
refined in 17 48,i9^ w h ere the unitarity-based decomposition of multi-loop integrands has been 
addressed as a polynomial decomposition problem, and systematized within the multivariate 
polynomial division algorithm. Accordingly, any generic multi-loop integral with n denominators, 
I12 ...n = J d d qi---d d q m I\2...n , with I\2... n = A r i 2 ...n/(T ) i • • • D n ), can undergo an integrand 
decomposition by means of successive polynomial divisions (modulo Grobner basis) between the 
numerator and the denominators, see Fig.2. The result of the decomposition reads as, 

J _ A 12 ...n A 2 ...n A 12 ... 71 —1 A n Ai 

12 " n ~ D\...D n D 2 ... D n " T>i... D n -1 " Dn " D~i 

where Aare the remainders of the iterated divisions (w.r.t. the Grobner basis of the ideal 
( Di ,... ,Dj )). Each residue Ais a polynomial in the components of the loop momenta not 
constrained by the cut D t = ... = Dj = 0. Therefore, by integrating both sides, one obtains the 
decomposition of the original integral Xi 2 ...n in terms of independent integrals. The integrand 
decomposition (1) implies that, exactly as it happens for the tree-level amplitudes, also the 
integrands of multi-loop amplitudes can be decomposed in terms of independent building blocks 
simply by partial fractioning ! 

While in the one-loop case the independent integrals are analytically known, in the multi¬ 
loop case, their classification and evaluation is an open problem. 


3 Differential Equations and Feynman Integrals 


The method of differential equations (DEs) 20,21,22 , reviewed in 23,24 ’ 25 ^ i s one of the most effective 
techniques for computing dimensionally regulated multi-loop integrals. 

In fact, any f-loop integral X is a homogeneous function of 
external momenta pi and masses m;, whose degree 7 = 7 (d,£) 
depends on the space-time dimensions d = 4 — 2e, on the number 
of loops £, and on the powers of denominators. Therefore, one 
can write the Euler scaling equation, 

( J2 Pi ' dr P' + d m{) X = ^) x > ( 2 ) 

* j 



Figure 4 - The three-loop ladder 

where d x = d/dx. Euler relation can be engineered to show box diagram, with one off-shell leg 
that Mis obey linear systems of first-order differential equa- ( daslictl linc ) 
tions (DEs) in the kinematic invariants, which can be used for 

the determination of their actual expression. By establishing an analogy between Schrodinger 
Equation in the interaction picture (in presence of an Hamiltonian with a linear perturbation) 
and systems of DEs for Feynman integrals (whose associated matrix is linear in e) 26 , we have 
recently proposed an algorithm to find the transformation matrix yielding to a canonical sys¬ 
tem 27 , where the dependence on the dimensional parameter e is factorized from the kinematic. 
In particular, we found that the canonical transformation can be obtained by means of Magnus 
exponential matrix 28 . The integration of canonical systems is simple, and the analytic proper¬ 
ties of its solution are manifestly inherited from the associated matrix, that becomes the kernel 
of the representation of the solutions in terms of repeated integrations. The latter in fact are 
the coefficients of a Magnus (or alternatively Dyson) series expansion in e. Magnus exponential 
is not unitary, as it happens in the quantum mechanical case, but the proposed method can be 
considered also inspired by unitarity. 



3.1 Applications 


We made use of Magnus theorem for the determination of non-trivial integrals, like the two-loop 
vertex diagrams for the electron form factors in QED and the 
two-loop box integrals for the 2 —> 2 massless scattering 26 , the 
two-loop corrections to the pp —»• Hj , as well as for evaluating 
the three-loop ladder diagrams for pp —> Hj (in the infinite top- 
mass approximation) 29 , see Fig.4. The latter is a formidable 
calculation involving the solution of a system of 85 Mis. In this 
case, after identifying a set of Mis obeying a linear system of 
differential equations in x = —s/m 2 H and y = —t/m 2 H , by means 
of a Magnus transform, the system can be brought in canonical 
form, reading as, 

dl(x, y) = e A(x,y) l(x,y) , 

where X is the vector of Mis, and df = d x fdx + d y fdy. The matrix A is purely logarithmic, 
A(x, y) = ai ln(x) + 02 ln(l - x) + as ln(y) + 04 ln(l - y) + <35 ln(x + y) + ae ln(l - (x + y)) , 
where the a* (* = 1,..., 6) are 85 x 85 matrices whose entries are just rational numbers. The 
logarithmic form of A trivializes the solution, which can be written as a Dyson series in e, where 
the coefficient of the series are combinations of Multiple Polylogarithms with uniform weight 
(where the weight increases as the order in e does). 

Boundary conditions are determined by 
imposing the regularity of the solutions in spe¬ 
cial kinematic configurations. Surprisingly, to 
fix the boundary values of all 85 Mis, only 2 
simple integrals have to be independently pro¬ 
vided. 

Also, we have been considering the mixed 
EW-QCD corrections to Drell-Yan production 
at NNLO, whose representative diagram is de¬ 
picted in Fig.5. Also in this case, Magnus ex¬ 
ponential can employed to reach a canonical 
system for the 48 Mis drawn in Fig. 6, in the 
variables x = —s/niy and y = — t/rriy , with 
V = W,Z 30 . 

4 Conclusions 

In this contribution, I have analyzed the algebraic patterns underlying the structure of scattering 
amplitudes in gauge theory. Unitarity plays a central role in the context of evaluating scattering 
amplitudes. It not only inspired methods to perform the amplitudes decomposition, by means 
of unitarity-cuts, but it also suggested a technique for the evaluation of master integrals, by 
means of matrix exponentials, similar to the unitary time-evolution in quantum mechanics. 



Figure 6 - Master integrals for qq —> E + £~ at two-loop. 
Plane lines stands for massless particle, while dashed 
lines stands for massive particles. 



Figure 5 - Representative two-loop 
box diagram for QCD-EW correc¬ 
tions to qq —> £ + l~. 
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